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Abstract 



We investigate the angular distributions of the rare B decay, B — > K*{^ 
KTT)i'^i^ , in general supersymmetric extensions of the standard model. We con- 
sider the new physics contributions from the operators 07.8,9,10 in small invariant 
mass region of lepton pair. We show that the azimuthal angle distribution of 
the decay can tell us the new physics effects clearly from the behavior of the 
distribution, even if new physics does not change the decay rate substantially 
from the standard model prediction. 
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I. INTRODUCTION 



Rare decays of B meson, e.g. b ^ and b — > si'^i'^, are the most suitable candidates 
for study of new physics beyond the standard model (SM). Since their branching fractions are 
usually very small within the SM predictions, they can sometimes show up new physics with 
unexpectedly large values of decay rates. The decay b —>■ s'j, which has been already measured 
by the CLEO Collaboration IQ , has shown that there is not much extra parameter space for its 
branching fraction from new physics . The decay b — * si~^i~ suffers from severe backgrounds 
of J/ip and ip' resonance contributions in the measurement of branching fraction, and so it 
may not be easy to uncover new physics cleanly . However, a number of methods have been 
discussed to detect new physics through the details of the decay process, such as differential 
distributions or polarization effects [Q-Q. 

Recently a new method utilizing the angular distribution of i? ^ K*{^ KiT)i~^i~ has been 
proposed |0 . Imagine the decay configuration when K* is emitted to the direction of +z and 7* 
is emitted to the opposite direction in the rest frame of B meson. Here 7* is off-shell photon and 
it further decays into i^i^ , and K* subsequently decays into Kn. If we ignore small mixture of 
the longitudinal component, the angular momentum of K* is either = +1 or = — 1, and 
the corresponding production amplitude is proportional to Cfji or Cji, respectively. Suppose 
the final K meson is emitted to the direction of {6k, <p) in the rest frame of K* ., where 9k is a 
polar angle and (f) is an azimuthal angle between the decay plane of {Kix) and the decay plane 
of In the low invariant nif+f- region, electromagnetic operator terms are dominated 

and the decay amplitude for the whole process is proportional to 

A C7iexp(-i0) + B C7/jexp(+z0) + C , 

where A, B and C are the real functions of the other angles. 

In this new method, we can distinguish the new physics contribution from that of the 
SM even if the branching fraction of the decay is similar to the prediction of the SM: In the 
B K*'y decays the probability of B meson decaying to left-handed (or right-handed) circular 
polarized K* is proportional to ICrxp (or ICy/jp), and therefore the polarization measurement 
of K* and 7 is useful for extracting the ratio of ICyil/ICy/jl. Even though the polarization of 
high energy real photon cannot be measured easily, we can still get some useful information 
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through the azimuthal angle distribution in the low invariant mass region of dileptons for 
B K*{~^ K7i)i~^i~ decay. This is because the decay products of K* and the virtual photon 
7* are responsible for this polarization measurements. In the SM, the operator O^l is dominant 
and the operator Ojr is suppressed by 0{ms/mi,). In this case, the angular distribution of the 
decay products is a flat function of the angle in the small lepton invariant mass region. If 
there is new physics contribution, the contribution of both operators can be equally important 
1^. We can distinguish the new physics signal easily from the angular distribution of the 
decay B K*{-^ K7[)i^i~ , while the measured branching fraction for B X^j can still be 
accommodated. 

In this paper, we extend this method to calculate the angular distribution of i? — *■ 
K*{Ktt)[^[~ in generalized supersymmetry models (gSUSYs). In addition to the operator 
OiR, we also consider the new operators O^r and Oiqr. In the following section, we derive the 
general formula for the various angular distributions including the new operators. Section III 
is devoted to the numerical analysis and discussions in gSUSYs. 



II. ANGULAR DISTRIBUTIONS OF THE DECAY B ^ K*{^ K + tt) + i+ + 

We start with the general effective Hamiltonian for the corresponding h sC^C.~ decay |^, 

n,ff{h s^+n = -^Vt,v:- + c,rO,r) . (i) 

The operators Oj relevant for us are 

OjL = ^{sLaMF'", (2) 

OrR=^{sRa,MF'^^, (3) 

OsL = ^ m,SLa^''T%R G^, (4) 

OsR = ^ m,SRa''''T%L G%, (5) 

0,L = {sb)L{U)v, (6) 

Om = {sb)R{ii)v, (7) 

OioL = {sb)L{ii)A, (8) 

OioR = {sb)R{ie)A, (9) 
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where in addition to the SM operators Ojl, Osl, Oql and Oiol, we include new operators Ojr, 
OsR, OgR and Oiqr. The new physics effects can originate from any of the above operators. 
The operator Ojr has aheady been introduced in Ref. [0, while the operators Oqr and Oioi? 
are newly introduced in this paper. The operator Og is a chromo-magnetic dipole operator. 
Relegating the details to Ref. [0, we introduce the helicity amplitudes for the decay 



B K*i 



which can be expressed by 



Ho' 



(ai + clVX 

P-L hX 

pi pi 
{clr + CRy/Xj , 

pi pi 



(10) 



with P = pk +Pn, L = p+ + p_, p = VP^ / = VL2, and A = (m|j - p- l^y/A - pH\ And 



Or, bu, cr and ai, &l, cl are given by 

aL = -C-j- \2{P-L)g+ + L^ 
bL = 



+ 9- 



+ 



{Cg- — Cio_)/^2 



2mb 



Cio+)g 



nib 



(11) 
(12) 

(13) 



Cr- 2{P-L)g+ + L\g+ + g-) 



2mh 



bR = -2Cr.{g+ - L'h) - 



Cr = -2C-j+g+ + 



(C9+ + Cw+)9 
nib 



nib 



(14) 
(15) 
(16) 



where the form factors g, (7+, g^, /, h and a"*" of 5 ^ K* decay are defined in Refs. ||7|. [Tl|JT2 
We introduced the Wilson coefficients C7_,7+,9_,9+,io-,io+ as 
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Ct- — Ctr — CrL, (^7+ = C7R + C7L, 

Clo- = CioR — CioL, Cio^ = CioR + C'loL- 

Using the above helicity amplitudes, the angular distribution of B 

expressed by, 



K*{^ K7i)i+£- is 



a 



* 12 

ts\ 



dp^ dp d COS dxd COS e+d(t> 64 x %{2T:fm%P[{p^ - mj^.)'^ + mj^.r^.] 

X [Acos^eKsin^e+(\H^\^ + \H^\^) 
+ sin^ 9Kil + cos^ e^)i\H^,\' + \H^,\' + \H^,\' + \H^,f) 
-2 sin^ Ok sin^ 6+ [cos 24)Re{H^^H^l + H^iH^l) 

- sin 20 Im{H^^H^l + H^^H^l) 

- sin 2^if sin 2^+ [cos (f)Re{H'^^H^* + i/:^ii/o^* + iZ+ii^o^* + H^l^H^*) 

- sin Im(H^^H^* - H\H^* + H^^H^* - H^iH^*) 



-2sin2 eKCOse+{\H 



R |2 
+1I 



1^1 



i? |2 



L |2 
+1I 



+ 1^-1 



L \2\ 



+2sin9+sin29K [cos 4)Re{H^^H^* - H^^H^* - H^iH^* + H^iH^*) 
- sin Im{H^iH^* + H\H^* - H^^H^* - H^^H^*)] } . 



(17) 



Here we introduced the various angles as: 6k is the polar angle of the K meson momentum in 
the rest system of the K* meson with respect to the hehcity axis, i.e. the outgoing direction 
of K*. Similarly 9+ is the polar angle of the positron in the 7* rest system with respect to 
the helicity axis of the 7*. And is the azimuthal angle between the planes of the two decays 
K* Kn and 7* ^ 

If we integrate out the angles 9k and 9+, we get the (p distribution 



dr 

d^ 



I 



a 



'GWK*K.VXp'ml\Vt,V*\' 



9 X 16(27r)8m|/2[(p2 _ mj^.y + m^.r|^ 
\H^\'^ + \H^i\'^ + \Hti\'^ - cos2(t)Re{H^^H^l + H^^H^l) 

+ sm2(t)Im{H^^H^l + H^^H^l)} dp^df . 



(18) 



Similarly, we can get the 9k and 6'+ angular distributions as following: 
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+ \H^\') + sin^ Ok {\H^,\' + \H^^\' + \H^^\' + \H!:,\')} dp'dl' , (19) 

and 

TjR\2 I \ttL\2\ 

dcose+ J 6x6A(27rfmlP[(p^-mj,.y + mlSl.]^ t^+U^o I + l^o N 

+ (1 + {\H^,\' + \H^,\') + (1 - cos^+)^ {\H^,\' + \H^,\')} dp'dl' . (20) 

Taking the narrow resonance limit of K* meson, i. e. using the equations 

- 487r ' 

rll"^o (p2_^2^,)2 + ^2^^rl,. = "'^^^ - "^-^ ' 

we can perform the integration over p^ and obtain the double differential branching fraction 
with respect to dilepton mass squared P and angle variables, 

dI^ = ^^ 384^^^Wl^'^*^^*^' 2^1'^°' 

+ + \Hl^\^ + - cos2(l)Re{H^^H^l + Hl^H^l) 

+ sm2(l)Im(H^^H^l + H^^^H^l)} , (22) 



rL\2\ 



+sin"fe(iffsr+iw?ii'+i//ir+ii/f:ir)} , (23) 



+ (1 + COS^+)2(|if^,|2 + + (1 - C0S^+)2(|i/f J2 ^ |if^j2)| ^ (24) 



and 



dP ~ ""sSAn^^' msH^ '''''''' I'^^o I - I-+1I 

+ |i/o''|' + |^^il' + |i^-il'} , (25) 

where tb is the life time of B meson, and we replace all p by ttik* due to the 6 function. 
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Finally, to eliminate the constant factors in Eq. 



we define the normalized distribution, 



s = 



dBr 
dPdcf) 

271 



I 



dBr 



dP 

cos2(j)Re{H^^H^* + H^^H^\) - sm2(f)Im{H^^H^l + H^^H^l) 



(26) 



where s = t^/m?^. The distribution r{(f),s) is the probability for finding K meson per unit 
radian region in the direction of azimuthal angle 0. Therefore, the normalized distribution 
r(0, s) oscillates around its average value given by 2^ ~ 0.16. 



III. GLUINO MEDIATED FLAVOR CHANGING NEUTRAL CURRENT AND ITS 

NUMERICAL ANALYSES 



In this section we consider the flavor changing neutral current (FCNC) in generalized su- 
persymmetric models (gSUSYs). In gSUSYs, the soft mass terms for sfermions can lead to 
potentially large FCNC 0. In the mass-insertion-approximation (MIA) [^, one chooses a 
basis for fermion and sfermion states, in which all the couplings of these particles to neutral 
gauginos are flavor diagonal. Flavor changes in the squark sector are provided by the non- 
diagonality of the sfermion propagators, which can be expressed in terms of the dimensionless 



parameters {6l 



IMN, 



iSlj)MN 



{''^ij)'MN 

m 



{M,N = L,R) 



(27) 



where {Tn1j)\.j]^ are the off-diagonal elements of the q mass squared matrix that mixes flavor 
i,j for both left- and right-handed scalars, and m is the average squark mass. The expressions 
for the Wilson coefficients at Mw scale due to the FCNC gluino exchange diagrams are 



w 



^SUSY//.^ 
'-'8L V^-^W 



9V2G^m2At 
__Tras__ 

V2GFm^Xt 



(^y,,M3(x) + (4),,;^M,(x) 



^^Msix) + 3M4(x)) + (52^3)^^ ^ (^Mi(x) + 3M2(x; 



'23 



LL 



Plix) 



{2i 



and 
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^SUSY,'/.^ 



V2GFm^X 



with Xt = V^lVtb- The functions Mi^3(x) and Pi(x) are defined as 

1 + 4x - + (4x + 2x^) Inx 



Mi(x) 

M2(x) = -x' 
MAx 



2 5 - 4x - + (2 + 4x) In X 
2(1 -x)4 ' 
-1 + 9x + Qx^ - 17x^ + (ISx^ + 6x^) In x 

12(x - If 
— 1 — 9x + 9x^ + x^ — 6(x + x^) Inx 



6(x - 1) 
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^ , , 1 - 6x + 18x2 - lOx^ - 3x^ + 12x3 In x 

-T = ; 77 , (30) 

^ ^ 18(x-l)5 ' ^ ' 

where x = rn^/fh^ and is gluino mass. 

In addition to the above gSUSYs contributions, the usual SM contributions Cj^, C'Il', C*!!^, 
and Cfi^ are aheady known for years, which we will not show here. Please look at Refs. for 
details of the SM contributions. Including the QCD corrections, we get the Wilson coefficients 
at rrib scale as 

CtlK) = -0.31 + o.erc^susY^^^) ^ o.ogcf^sY^^^^ ^ (3^^ 

Cr^K) = 0.67CSUSY^^^) ^ 0.09CSUSY(M^) , (32) 
C9L(m,) = Cff (m,) + CfusY(M^) , (33) 
C9K(mfe) = CfU"^(Miy) . (34) 

Here, —0.31 in Eq. (|3T|) is the SM value of C-jLijnb)- New physics contributions Cj^^'^(Mvy), 
Cf^^'^(Mvi/) come from Eqs. ( PB| , pUD . For Cio, there is no new physics contribution. 



CioL(m,) = Cfo'^(mb) , (35) 
Cioi?(m,) = . (36) 

Then we can get the complete formula for angular distributions of i? — > K*i^i^ , Eq. 
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The operators 0-jl and O-jr contribute to the rare radiative decay h — > 37. Their Wilson 
coefficients have been constrained by the experimental measurements of the decay. The decay 
width for inclusive 6 — > 57 decay is given in terms of the operators 0-jl and 0-jr. It is convenient 
to normalize this radiative partial width to the semileptonic decay h cev in terms of the ratio 



r(6- 


si) 


_ 6 \VtlVtb\'^ aem 


ICjLirrib) 




r(6- 


^ cev) 






{mb)g{mjmh) 



where the functions f{x) and g{x) are phase space and QCD correction factors |jT3[, respectively. 
The b s'-f branching fraction is obtained by 

Bn{b ^ 57) ~ Bn{B XJz/)exp. X i? ~ (0.105) X R. (38) 

For ETZ{h —>■ 57), we use the present experimental value of the branching fraction for the 
inclusive B — > Xsi decay, 

Bn{B X,7) = (3.15 ± 0.35 ± 0.32 ± 0.26) x 10"^ (39) 

Constrained by this experiment, we derive from Eq. (|37D 

\C7L{mb)\^ + \Cm{mt)\'' = 0.081 ± 0.014. (40) 



In the numerical calculations, we use the form factors calculated in Ref. |jTT[. They are listed 



in Table |I| for zero momentum transfer. The evolution formula for these form factors is 

where l"^ = {pi+ +p^-)^. The corresponding values ai and a2 for each form factors are also 
listed in Table 1. 

The decay B K* + Jji\}{—^ t^i~) is a possible background for our B K*l^l~ decay 
at the J I'll) resonance region, so as the ip' , etc. Therefore, only the low invariant mass region of 
the lepton pair is good for clean measurements. The helicity amplitudes are dominated by the 
two coefficients C^l and Cju in the region of low invariant mass, as given by 

H^'^ ~ . (42) 



In the small invariant mass limit s <ti 1, r(0, s), defined in Eq. (^), is approximately written 

as, 

, « 1) ^ _ 1 1 + ,0,2^ ^ ^- ^^ _ I ^y^^i ^ . (43) 



In the SM case, Cyj? ~ and therefore the above approximate formula is reduced to 

r(0, s ^ 0)sM ^ (44) 

In Fig. 1 we can see that it is almost a constant distribution of in the small s region. As 
s increases, the contributions from the operator Og and Oio makes (~ —cos 20) behavior. 
However, the new physics contributions can give quite different distributions depending on the 
model, and we can probe new physics efficiently. Here we discuss the gSUSYs contribution to 
the distribution r(0,s). For simplicity, we assume that \5lr\ = \5rl\ and = \Sb.r\, and 
consider two cases: LR mixing dominating case and LL mixing dominating case. 

First we consider the LR mixing dominating situation, i.e. ~ \Srr\ ^ \^lr\ ~ \Srl\- 
Fig. 1^ shows the distribution r(0,s) for 6lr = Srl = \Xt\ case with x = 0.3, m = 960 GeV. 
This corresponds to Cjl = 0.017 and Cjr = 0.333. Since Cjl Cjr and both are real, the 
approximate formula (^3]) becomes 

r(0, s<l)~^(l + cos20 ^1. (45) 

This (~ cos 20) behavior is shown in Fig. ^. On the other hand. Fig. ^ shows the distribution 
r(0, s) for 6lr = —Srl = \Xt \ case with x = 0.3, m = 960 GeV. This corresponds to C^l = 0.017 
and C-m = —0.333. In this case the approximate formula (^) becomes 



s < 1) ~ <i 1 - cos: 



(46) 



CrR 

This (~ — cos 20) behavior is shown in Fig. ^explicitly for the small s region. In Figs, ^and 
we used the same values of x and m: The branching fractionss of both b ^ s'y and b — si~^i~ 
are unchanged for the two situations, and we cannot separate these two situations by using only 
branching fraction measurements. However, we can see that the angular distributions, shown 
in Figs. 2 and 3, can easily distinguish the relative sign of Slji and 6jiL- 

For the LL mixing dominating case, i.e. \Slr\ ~ \Srl\ ^ ~ \Srr\, we also show two 
cases. First we choose 6ll = Srr = e*5, with x = 0.8, rh = 250 GeV. This corresponds to 
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CjL = —0.16 + i 0.15 and Cjr = +0.15 + i 0.15, i.e. Cjr/Cjl ~ case. Using this set of 
parameters, the formula (|43|) becomes 



Fig. H shows this (~ sin 20) behavior clearly. On the other hand. Fig. 5 shows the distribution 
r(0, s) for 5ll = —Srr = e*? case with x = 0.8, fh = 250 GeV. This corresponds to Cjl = 
—0.16 + i 0.15 and Cjr = —0.15 — i 0.15, i.e, C-jr/C-jl ~ case. The approximate formula 
becomes 



The (~ — sin 20) behavior is shown in Fig. 5. From Figs. 4 and 5, we can see that the angular 
distribution can distinguish the relative phase between 5ll and 5rr easily, even if we use the 
same values of parameters, x and ffi. 

The polar angle distribution functions in Eqs. (PB| , |^ ) depend only on the modular square 
terms of the helicity amplitudes, which give the decay width of the semileptonic decay. If the 
branching fraction is fixed by experiments, these two angle distributions cannot distinguish 
new physics contribution from the SM. On the other hand, they can serve as a double check of 
whether the branching fraction is different from the SM predictions. 

In conclusion, we have calculated the angular distribution of the rare decay, B — > K* + (— *■ 
Kti) + + l~ in general supersymmetric extensions of the standard model. The azimuthal 
angle (0) distribution in gSUSYs can be quite different from that of the SM, while the measured 
branching fraction for B Xs'y can be accommodated within the standard model prediction. 
In the standard model it is found to be almost a constant under the variation of the angle 
in small invariant mass region, while in gSUSYs the distribution can show (~ ± cos 20) 
or (± sin 20) behavior depending on the gSUSYs parameters. We showed that the angular 
distribution of the decay can tell us the new physics effects clearly from the behavior of the 
distribution, even if new physics does not change the decay rate substantially: We would be 
able to tell the relative phase between the mixing parameters 6lr and 6rl (or 6ll and Srr), 
even though the decay rate of gSUSYs were exactly the same as that of the SM. 




(47) 




(48) 
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TABLES 



TABLE 1. Form factors in zero momentum transfer and parameters of evolution formula [11]. 





9 


/ 


«+ 


a_ 


5+ 


9- 


h 




0.063 


2.01 


-0.0454 


0.053 


-0.3540 


0.313 


-0.0028 


o\ 


0.0523 


0.0212 


0.039 


0.044 


0.0523 


0.053 


0.0657 


02 


0.00066 


0.00009 


0.00004 


0.00023 


0.0007 


0.00067 
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FIGURES 




FIG. 1. The distribution r((p, s) for the SM case, where C7R ~ 0. 





FIG. 3. The distribution r{(p,s) for S^r = —6rl = \Xt\ case, with x = 0.3, rh = 960 GeV. This 
corresponds to Cjl = 0.017, Cjr = —0.333. 



r 




FIG. 4. The distribution r{(p,s) for 6ll = 6fi.r = e*4 case, with x = 0.8, rh = 250 GeV. In this 
case, Cyr/C-tl ^ e'^'^. 
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